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In this work we have studied the effect of target thickness on relativistic self-induced transparency (RSIT) 
and found out that for subwavelength targets the corresponding threshold target density (beyond which target 
is opaque to incident laser pulse of a given intensity) increases. The accelerating longitudinal electrostatic field 
created by RSIT from subwavelength target is used to accelerate the ion bunch from a thin, low density layer 
behind the main target to ~ 100 MeV. A suitable scaling law for optimum laser and target conditions is also 
deduced. The word ‘via’ in the title signifies the fact that we are interested in acceleration of ions from the layer 
placed behind the target. It is also being observed that as per as energy spectrum is concerned; an extra low 
density layer is advantageous than relying on target ions alone. 


I. INTRODUCTION 

The laser based particle acceleration have been shown a 
great deal of research interest since last couple of decades and 
still continues to do so. The researchers in this field around 
the globe are working towards the understanding and the de¬ 
velopment of the efficient acceleration mechanism both on ex¬ 
perimental and theoretical front. The development of efficient 
particle beam source promises remarkable applications in var¬ 
ious areas of applied sciences [1-6]. 

The laser and target conditions play paramount role in de¬ 
ciding which acceleration mechanism would prevail. Among 
all the acceleration mechanisms the Target Normal Sheath Ac¬ 
celeration (TNSA) is studied extensively both on experimen¬ 
tal and theoretical front [7-13]. The weak scaling of proton 
energy with laser intensity (~ s/Tq) is a severe limitation of 
the TNSA as per as GeV protons are concerned [7]. When 
the target thickness is comparable to skin depth of the plasma, 
then after TNSA another mechanism namely Break-Out Af¬ 
terburner (BOA) takes over. In this stage, the hot electron 
plasma expands sufficiently by laser pondermotive force so 
that target becomes underdense. The laser penetrates the tar¬ 
get and strong electric field is induced through Buneman insta¬ 
bility, this growth of instability results in conversion of elec¬ 
tron energy to ion energy and thus, ions are accelerated to very 
high energies [14-16]. 

When the laser intensity increases to ~ lO^** W/cm^ then 
radiation pressure acceleration (RPA) starts to dominate the 
acceleration process [17-24]. Depending on the target thick¬ 
ness two mechanisms can occur under RPA viz. Hole Boring 
(HB) for thick targets and Light Sail (LS) for thin targets. The 
RPA is mainly governed by circularly polarized laser pulse so 
that j X B heating can be avoided and electrons are pushed 
deep into the targets. In HB radiation pressure drives material 
ahead of it as a piston but does not interact with the rear sur¬ 
face of the target [25, 26]. However, in alternate scenario if the 
target is sufficiently thin for the laser pulse to punch through 
the target and accelerate a slab of plasma as a single object is 
referred as LS regime [27, 28]. A wonderful review on this 
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contemporary field of laser induced ion acceleration is pre¬ 
sented by Macchi et. al. [29] in which all the ion acceleration 
mechanisms are covered in depth. 

In order to efficiently accelerate the ions to ultra-relativistic 
energies the use of multilayer or multi-species target is also 
examined by various researchers [20, 30-33]. In Ref [34] 
authors have shown the acceleration of protons to ~ 210 MeV 
using double layered thin foil by the laser with peak inten¬ 
sity 2.7 X 10^^ W/cm^. In this so called Directed Coulomb 
Explosion (DCE), the high intense laser beam interacts with 
very high density Al*^+ plasma (~ 400 nf) which expels all 
the electrons from the target leaving behind ions, which un¬ 
dergoes Coulomb explosion, the adjacent proton layer of 30ne 
then gets accelerated by the electrostatic fields created by the 
DCE [20]. 

Each mechanism defers from another on the basis of laser 
and target conditions used. The relativistic self induced trans¬ 
parency (RSIT) can also be perceived as one of the mechanism 
for acceleration of ions, this exploits the fact that the thresh¬ 
old plasma density for laser penetration, increases as a con¬ 
sequence of relativistic mass effect of the plasma electrons. 
RSIT was initially reported in pioneer work of Refs. [35, 36], 
since then it has also drawn considerable interest in ion accel¬ 
eration fraternity. RSIT can prevent the RPA in thin targets, 
and hole-boring in thick targets. Eurthermore, RSIT can re¬ 
sults in hot electrons and thus allowing for more energetic ion 
acceleration [37-39]. Recently laser-driven Neutron source 
based on the relativistic transparency of solids is also being 
reported [40]. 

As per as RSIT is concerned, many of the theoretical dis¬ 
cussions seek the steady state solution in the purview of rela¬ 
tivistic cold fluid model. However, in this work we attempted 
to solve the set of dynamical equations in space and time 
which governs the response of semi-infinite plasma to in¬ 
tense laser fields. In order to study the RSIT for thin targets, 
we relied on relativistic stationary plasma model along with 
PIC simulations. In this work we would be using RSIT phe¬ 
nomenon to accelerate the protons to very high energies using 
a secondary thin and low density layer behind the main target. 
The protons from this secondary layer will then be accelerated 
as a bunch in the electrostatic field created by RSIT mecha¬ 
nism in main target. The dependence on the target thickness 
and density is also being studied and it has been observed that 
the threshold density for RSIT increases for subwavelength 
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targets. The objective of this article is to understand the for¬ 
mation of electrostatic fields by RSIT mechanism and its use¬ 
fulness to accelerate protons. In section II we briefly discuss 
the RSIT mechanism along with a theoretical and simulation 
model for the same. The detailed formulation of the dynam¬ 
ical equations are discussed in Appendix A. Results and dis¬ 
cussions are presented in section III followed by the conclud¬ 
ing remarks in section IV. 


II. THEORY AND SIMULATION MODEL 


The propagation of electromagnetic (EM) waves in plasmas 
have been studied quite extensively and vast literature on the 
topic can be found [41, 42]. The dispersion relation for EM 
waves propagation in plasmas is expressed as 

o? = (Op+k^c^, ( 1 ) 

where O), and k are respectively the frequency and wave prop¬ 
agation vector of EM waves, c is speed of light in vacuum 
and (Op = \JIEaine is the natural frequency of the plasma 
oscillations. Here e, and are respectively the electron 
charge, electron mass and electron density. 

It can be understood from Eq. 1 that EM waves can not 
propagate beyond the point where (Op > (0 because it meant 
the wave vector to be imaginary, which is physically equiva¬ 
lent of attenuation of EM waves. In terms of the plasma den¬ 
sity it means that EM wave can not propagate in plasma be¬ 
yond a critical density ric which is defined as iic = EanieOp'/e^ 
(it is the density for which the plasma frequency matches to 
that of EM wave), however for densities the EM wave 

can propagate without much of attenuation. 

The interaction of very intense laser fields can change the 
above mentioned criteria about the critical density. In such 
situations the relativistic mass of the electrons need to be 
taken into account, which in principle increases the thresh¬ 
old density (n^) for the propagation of laser for a given inten¬ 
sity. This effect is refereed as relativistic self induced trans¬ 
parency (RSIT). The interaction of ultra intense laser pulses 
with plasma can be modeled by solving the wave propagation 
equation in Couloumb gauge, continuity equations, posisson’s 
equation along with relativistic Lorentz force equations. The 
analytical treatment of RSIT based on relativistic stationary 
plasma model (cold fluid theory with steady-state solutions) 
is being vastly discussed in the litretaure [43-45]. However, 
in this work we have attempted to solve the set of dynam¬ 
ical equations in space and time which governs the plasma 
response to intense laser fields. We would restrict ourselves 
to the one dimensional scenario where laser is considered to 
be propagating along z direction, in view of this the spatial 
derivatives along z directions is only considered. The detailed 
derivation describing the plasma response (electron + ions) to 
intense laser fields is discussed in Appendix A. In the absence 
of the ionic motion the closed set of equations describing the 
plasma response reduces to (Appendix A), 
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Ye = (p|)2 (8) 

where, Z is atomic number of the ionic species, rig and n, are 
electron and ion densities in units of critical density, p® is lon¬ 
gitudinal electron momentum, is longitudinal electrostatic 
held created by charge separation, jg is relativistic factor of 
electron, and a is vector potential associated with laser pulse. 
It should be noted that the motion of the electrons (Eq. 6) 
are governed by the electrostatic held and the ponderomotive 
force (djgldz) of the laser pulse (see Appendix A), which in 
turn manifests the lowering of the target density and hence as¬ 
sisting the laser propagation deeper into the target, as can be 
inferred from Eq. 5. Eurthermore, it can also be seen from Eq. 
5 that in rig/jg —> 0 then it will be having a solution of EM 
wave propagation in vacuum. In order to have a RSIT regime 
the condition should be rigljg < 1 or the laser would be re¬ 
flected from the overdense plasma. As can be seen from Eq. 8 
that one in principle needs to solve for p\ to have the knowl¬ 
edge of the jg and hence the threshold plasma density where 
RSIT ceases to exist. In this context one can infer the fact the 
electron heating (p^ plays quite a crucial role in RSIT. 

In order to study the interaction dynamics in RSIT regime, 
we have solved these equations numerically in space and time 
with following initial/boundary conditions, 

5cos(f)ex+ \/l — 5^sin(f)ey 

( 9 ) 

where, T is EWHM pulse duration, and 5 is a constant which 
is either 1 for linearly polarized (LP) laser pulse or 1 /v/2 for 
circularly polarized (CP) laser pulses. In this study we will 
be using only CP laser pulses. In order to benchmark the re¬ 
sults of this analytical framework, we have simulated the in¬ 
teraction of 5 cycle Gaussian laser pulse (ao = 20) with semi- 
infinite plasma. The laser is allowed to incident normally from 
vacuum (z < 12A) to a semi-infinite plasma slab (z > 12A) 
with density n,- — ng = Irig. 

In Pig. 1 we have compared the longitudinal electrostatic 
held (Ej.) as calculated by numerically solving Eq. 2 to 8 
with ID PIC simulations, and as can be seen from the same 
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FIG. 1. (color online) The component of the laser field along x direction (£^) and electrostatic field (fij) created by charge separation is 
presented by numerically solving Eqs. 2-8 (left panel) and compared with ID PIC simulations (right panel). 


that the outcome of theoretical formulation is reasonably in 
good agreement with PIC simulations. However, in order to 
study the laser interaction with target of finite thickness, one 
needs to incorporate the plasma expansion in vacuum [46] 
with above set of equations, which is currently beyond the 
scope of the present work, and we would like to reserve the 
same for the future. 

Nevertheless, in order to have some qualitative understand¬ 
ing of the effect of finite target thickness and corresponding 
threshold density (beyond which the target will not be trans¬ 
parent for a given laser pulse profile), one can in-principle 
seek steady state solution in the transparency regime (other 
case might be the reflection of incident laser pulse by over- 
dense plasma). Except where otherwise stated, we will be 
studying the interaction of flat top laser pulse (aq) of duration 
6 cycles with rise and fall of 1 cycle each, with a target of 
thickness d and density n^. The pulse profile is considered to 
be flat top, such that the condition for RSIT (ng/yg < 1) can 
be met from the instant laser hits the target. Recently in Ref. 
[47], the relativistic cold fluid model with stationary solutions 
have been used for semi-infinite plasma and an expression for 
threshold plasma density has been reported. The threshold 
density for RSIT for high laser intensities is then given by 
[47], 


nth 
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We are mainly interested to have the relativistic transparency 
in a target with some finite thickness, such that the laser 
should pass through the target without much of attenuation. 
Though, the dependence of the target thickness is not ex- 



FIG. 2. (color online) Variation of threshold electron density (n,;,) 
with laser amplitude is presented. Curve ‘A’ shows the expression 
given by Eq. 10, points denote the results obtained by the ID fully 
relativistic PIC simulations for a target thickness of 1A (square), 0.5A 
(circle) and 0.25A (triangle). Curves ‘B’ and ‘C’ respectively show 
the qualitative estimate of threshold density for d = 0.5A and d = 
0.25A by using Eq. 11. 


plicitly appearing in Eq. 10, but we have tested this expression 
against a target thickness of lA using ID PIC simulation and 
results are presented in Eig. 2. It can be seen from Eig. 2, that 
the threshold plasma density for a target of thickness lA is in 
good agreement with Eq. 10, which in a sense indicates that 
the laser of amplitude qq would be sufficient enough to sweep 
out all the electrons from a target of thickness lA and having 
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maximum density of n,/, as predicted by Eq. 10. Total number 
of electrons present in a plasma slab (ID scenario) of density 
rig and of thickness d should be rigd, so the agreement of Eq. 
10 and PIC simulation will dictate that laser of amplitude qq 
would be sufficient to sweep out rigd electrons. We can extend 
this idea to subwavelength targets, and alternatively can write 
Eq. 10 such that for d =\X the threshold density is given by 
Eq. 10. The corrected expression for threshold density for a 
target of thickness d should read as, 


nth ^ 
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The results obtained by this simple scaling law given by Eq. 
11 is also presented in Eig. 2, and compared with PIC simu¬ 
lations for the same. As can be seen that the threshold density 
as predicted by Eq. 11 is reasonably in good agreement with 
PIC simulations. It should be noted that, in Pig. 2 the region 
below the data points for particular target thickness would be 
transparent and the region above the same would be opaque. 


III. RESULTS AND DISCUSSIONS 

In this paper we have used the 1D3V fully relativistic PIC 
code LPIC-h- [48] to carry out the studies. The dimension¬ 
less electric fields are normalized as ao = eE jmedrc, where 
CO is the frequency of the laser pulse, and E is electric field 
amplitude in SI units. Time and space are normalized by one 
laser cycle and wavelength respectively. We have modified the 
code to be able to study the laser interaction with multilayer 
target structure. The laser pulse is considered to be a flat top of 
duration 6 cycles with rise and fall of 1 cycle each, propagat¬ 
ing along the z direction and incidents normally on the target. 
Throughout the paper we will be using a simulation domain of 
length 60A (unless otherwise stated) and the Hydrogen plasma 
is considered to be located in region 15A <z<d, where d is 
the target thickness. Rest of the space is considered to be vac¬ 
uum. The target density will be chosen such that RSIT will 
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FIG. 4. (color online) Space-time dependence of electron density 
for target thickness 0.5A (a) and 0.25A (b) along with longitudinal 
electrostatic field (E^) for these two cases respectively are presented 
in (c) and (d). 


enable the laser to pass through. An extra layer of thickness 
0.2A and having a density of OArig is placed behind the main 
target. The density of second layer is kept low so that the laser 
can also pass through this layer as well, and the ions from this 
layer will then see the electrostatic field created by RSIT and 
will accelerate as a bunch. We have used the term ‘ions’ and 
‘protons’ interchangeably because we are only dealing with 
Hydrogen plasma in this paper. 

Before we start to study the effect of target thickness and 
density on RSIT and hence the accelerated ions, first we will 
see the difference between RSIT and LS regime of ion accel¬ 
eration. 
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FIG. 3. (color online) The spatial profile of electrostatic field (E^ and 
laser field (£^) is plotted at 26 t for three different target thickness 
having density Ing. 


A. How RSIT is different from LS regime of acceleration ? 

As per as LS regime of ion acceleration is concerned, it 
has been reported that the thickness of the target foil plays 
quite an essential role in determining the energetics of ac¬ 
celerated ions [27, 28]. If the ratio then the ac¬ 

celeration mechanism will be in LS regime and on the other 
hand the RSIT mechanism starts to prevail in the regime 
where the ratio ^ 1- Here, ao is the laser amplitude 
of circularly polarized laser and the parameter ^ is defined 
as ^ = K{ne/nc){d/X) with rig, ng are respectively the target 
density and critical density corresponding to laser wavelength 
A and d is the target thickness [27]. 


B. Formation of Electrostatic Field 

In the previous section we have discussed the basic theory 
behind the RSIT and elucidated the fact that the balance be- 
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FIG. 5. (color online) Velocity spectrum (a) of target ions and ions from secondary layer is measured at 90 t for the case qq = 20, d = 0.25A, 
tie = 3«c, 5nc and Iric- The energy spectrum (b) for He = Stic at different time steps is also presented. 


tween the ponderomotive force by the laser and the electro¬ 
static force by the charge separation is key and minimum cri¬ 
teria for RSIT mechanism. This balance of the forces plays 
paramount role especially when we have the plasma of the fi¬ 
nite length. The electrostatic force originated by the charge 
separation is very much depend on the thickness as well as 
on density of the target. As an example, in Fig. 2 we have 
presented the threshold density for RSIT as obtained by the 
PIC simulation for three different target thicknesses. It is be¬ 
ing observed that the thin target supports comparatively larger 
target density as compare to the thicker target as per as RSIT 
is concerned. 

In order to see the effect of the target thickness on the for¬ 
mation of the electrostatic field, we have simulated the laser 
interaction with a target of thicknesses 0.25A, Q.5X and IX. 
The laser conditions are same as in Fig. 2 (circularly polar¬ 
ized, 6 cycle flat top with rise and fall of 1 cycle) with gq = 20 
and target density for all the three cases is considered to be 
Trie- Figure 3 shows the spatial profile of the electrostatic field 
(a) and laser field (b) at 26t. Here, T is one laser cycle, and 
in vacuum laser propagates a distance of lA in time It. As 
can be seen from Fig. 3(a) that for the case of 0.25A the elec¬ 
trostatic field is almost flat in the region just behind the tar¬ 
get (15A <z<d), while for the target with thickness lA, the 
strength of the electrostatic field is not that prominent. This 
can be explained on the basis of the threshold density for dif¬ 
ferent target thickness as shown in Fig. 2. It can be read from 
Fig. 2 that for a target which is lA thick, the threshold density 
for ao = 20 lies slightly below the density used in this case 
(Fig. 3), however for 0.25A and 0.5A, it is higher than Iric- 
This bring the target with thickness 0.25A and 0.5A in RSIT 
regime, which in principle is responsible for sweeping all the 


electrons and creating strong electrostatic field. The trans¬ 
parency of the target can also be seen from Fig. 3(b), where 
the reflection of laser pulse for 1A thickness and transmission 
for 0.25A and 0.5A thick target is clearly seen. It should be 
also noted that the parameter ^ as defined earlier is about 5.5 
and 11 for respectively 0.25A and 0.5A target thicknesses, and 
hence the interaction with aq = 20 brings it in RSIT regime, 
however, lA thick target (i* ~ 22) is opaque for the laser with 
flo = 20. 

In Fig. 3, we have only shown the spatial profile of and 
Ex only at 26t, however in order to see the spatio-temporal 
evolution of the electrostatic field, the electron density and 
corresponding electrostatic field for the case with d — 0.25A 
and d = 0.5A are presented in Fig. 4. As we are using the 
same density (7nc) for thicknesses 0.25A and 0.5A, hence 
its expected that the electrostatic force would be weaker (in 
comparison with ponderomotive force) in case of 0.25A as 
compare to 0.5A, because effectively the target with thickness 
0.25A would be having less charge particles as compare to the 
case when thickness is 0.5A. The weaker electrostatic force in 
case of 0.25A manifests the stronger push on the electrons by 
CP laser pulse as can be seen in Fig. 4(b). On the other hand 
in case of 0.5A the ponderomotive force of the laser is not 
strong enough to sweep the electrons to large distances as can 
be seen from Fig. 4(a). However, in both the cases the elec¬ 
trons are pushed away from the target, lowering the density of 
the target and hence enabling the target to be transparent for 
laser to pass through. When the laser is passed through the 
target then the electrons are pulled toward the target because 
of the electrostatic force created by charge separation of elec¬ 
trons and ions. As can be seen from Fig. 4(a), the electrons 
oscillate at higher frequency as compare to their low thick- 
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ness counterpart as in Fig. 4(b). This can be understood from 
the lower excursion of the electrons in 0.5A case as compare 
to that of 0.25A case. It can be summarized as follows, for 
fixed density, thin (thick) target, will be having less (more) 
electrons, which results in weak (strong) electrostatic force to 
counter, and hence large (small) excursions of the electrons 
by ponderomotive force will take place, which in principle re¬ 
sults in low (high) frequency electron oscillations by electro¬ 
static force created by charge separation. The low frequency 
oscillations of the electron in case of thin target manifests the 
kind of constant electrostatic field which can be harnessed to 
accelerate the ions to very high energies. 


C. Acceleration 

So far we have discussed the formation of the electrostatic 
field by RSIT mechanism. Now we can exploit these fields 
to accelerate the protons to very high energies. Though it is 
clear that once the electrons are swept out from the target then 
because of the space-charge effect of the ions from the target 
will start moving in either direction, which would not yield 
in mono-energetic proton beam. In order to achieve the effi¬ 
cient acceleration by the electrostatic field created by RSIT, 
we have introduced a secondary layer of thickness 0.2A with 
density Q.lric just behind the main target (laser incidents on 
the main target first). The density of this extra layer is consid¬ 
ered to be low enough such that electrostatic fields created by 
RSIT mechanism is not affected by the presence of the sec¬ 
ondary layer. Secondly, as a consequence of low target den¬ 
sity, the laser will not have any direct effect on the energetics 
of the protons in this layer. 

The velocity spectrum of the ions from the main target and 
the ions from secondary layer for the case aq = 20, d — 0.25A, 
and He = 3nc, Stic and Trie as evaluated at 90 t are presented in 
Fig. 5(a). As can be seen from this figure that ions from sec¬ 
ondary layer are getting accelerated as a bunch, on the other 



FIG. 6. (color online) Variation of maximum energy of accelerated 
ions (from secondary layer) with density and thickness of primary 
layer for different laser amplitudes is presented. 


hand the ions from main target expands in either direction 
because of the space charge effect, resulting in large energy 
spread. This clearly shows the advantage of using an sec¬ 
ondary layer behind the target to have a mono-energetic ion 
bunch. In this case the target ions are only responsible for 
creating an electrostatic field which is exploited by an addi¬ 
tional thin layer of low density resulting in mono-energetic 
ion bunch. As per as the quality of the bunch is concerned, it 
can be seen that the spectra for the case Stic is much sharper 
than the one for the cases 3nc and Trie, which in fact signals 
toward the fact that there should be an optimum target density 
for a particular target thickness in order to have very efficient 
proton acceleration. The details regarding the optimum con¬ 
ditions for a given laser and target conditions is discussed in 
later part of this section. 

Furthermore, the energy spectrum of the accelerated pro¬ 
tons from the secondary layer for the case Wg = Stic as evalu¬ 
ated at different times are presented in Fig. 5(b). As can be 
seen that at early stage of acceleration, the protons gain energy 
at faster rate and later the rate of acceleration decreases, this 
is mainly because of the weakening of the electrostatic force 
over a period of time. It is being observed that the ions from 
secondary layer are being accelerated to ~ 100 MeV. 


D. Optimization 


As it has been pointed out in Fig. 5 that depending on laser 
and target conditions, there should be some optimum condi¬ 
tions in order to accelerate the ions most efficiently. In view 
of this we have presented the variation of maximum kinetic 
energy of the secondary ions with density and thickness of 
primary layer for different laser amplitudes in Fig. 6. The 
maximum value of the energy is evaluated at 90 t. It can be 
observed that there is an optimum density at which the ions 
can be accelerated most efficiently. This can be explained on 
the basis of the fact that for low densities overall charges itself 
is low enough resulting in weaker electrostatic field by pri¬ 
mary target. On the other hand, for higher densities the RSIT 
mechanism ceases to exist because of reflection of the laser 
pulse and hence resulting in poor longitudinal electrostatic be¬ 
hind the target. In view of this it can be understood that there 
should be some optimum target density for particular target 
thickness at which the interplay between laser ponderomotive 
force, and the electrostatic forces by charge separation, results 
an ambient condition for the ions from secondary layer to be 
accelerated. The optimization of the target thickness for RSIT 
is also being reported in an experimental paper by Henig et. 
al. [49]. 

The optimum density corresponding to particular laser and 
target conditions can be obtained from Fig. 6, e.g. for ao = 20, 
d — 0.25A the optimum density is about ~ 6nc- In Fig- 7 
we have presented the optimum density as calculated by PIC 
simulations. In this figure we have presented the optimum 
density corresponding to particular thickness and density of 
primary layer, and laser amplitude. It has been observed that 
the optimum density for particular laser amplitude (ao) and 
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FIG. 7. (color online) Optimum density of primary target is presented 
as a function of laser amplitude and its thickness. The solid line 
represents the scaling law for the optimum density as a function of 
laser amplitude and target thickness (Eq. 12). At these optimum 
densities the ions from secondary layer attains maximum energy. 


are accelerated most efficiently to few hundreds of MeV by 
electrostatic fields formed by RSIT. The role of target ions in 
this process is to create the electrostatic fields by RSIT mech¬ 
anism, which will eventually accelerate the protons from the 
thin and low density layer placed behind the target. A suitable 
scaling law connecting the optimum target density for maxi¬ 
mum energy, target thickness and laser amplitude is also being 
deduced. 

Though the results presented in this paper shed some light 
on the physics aspects of the acceleration of ions from sec¬ 
ondary layer in RSIT enabled accelerating fields in subwave¬ 
length target. However, the detailed quantitative analysis and 
optimization of the physical parameters is still warranted, we 
have reserved that for our future work. We also would like 
to extend the presented theoretical framework for the targets 
with finite thickness. It should be also noted that we have used 
ID PIC simulation to carry out these studies, and hence in or¬ 
der to experimentally realize these ideas one need to carry out 
full 3D PIC simulation of the same. 


target thickness (d) can be fitted with following scaling law, 
ne(ao,<i) = (0.07ao + 0.3)/(i (12) 

As can be seen that Eq. (12) fits very well with the predic¬ 
tions of the PIC simulations. The linear relation again hints 
towards the balance between the ponderomotive force ag) 
and the electrostatic force (°^ ngd) for maximum energy which 
in principle depends on the electrostatic force created by the 
charge separation in primary layer. In this work we have kept 
the thickness and density of secondary layer as 0.2A, and 0. Iwe 
respectively, in principle one should also have some optimum 
conditions on secondary layer, but we would like to reserve it 
for the future. 


IV. CONCLUDING REMARK S 

We have studied the RSIT from theoretical perspective by 
numerically solving the closed set of equations describing the 
plasma response to intense laser fields, and results are com¬ 
pared against the ID PIC simulations. The theoretical pre¬ 
scription tends to very well predict the strength of the longi¬ 
tudinal electrostatic field created by charge separation from 
laser ponderomotive force in semi-infinite plasma. The ex¬ 
isting theoretical framework need to incorporate the plasma 
expansion in vacuum in order to study the interaction of in¬ 
tense laser with target having finite thickness. In order to 
study the RSIT in finite thickness target, we relied on the sta¬ 
tionary plasma model to predict the threshold target density 
for a particular laser amplitude and target thickness and found 
that for subwavelength targets the threshold density for RSIT 
is increased. 

We have also discussed the role of target density and tar¬ 
get thickness on accelerating fields and it is being observed 
that the protons from the secondary layer behind the target 
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Appendix A 

1. In coulomb gauge the Maxwell’s equations will take the 
following form for vector potential A and scalar poten¬ 
tial (j), 


W^(j) = -p/eo 

(Al) 

1 d^A 1 fd^\ 

c2 dfi ~ 

(A2) 

J = —engVe + ZenjVi 

(A3) 

Ve = v"z-FV^ 

(A4) 

Vi = v‘z + vY 

(A5) 


where, all the symbols have their usual meaning. The 
index e and i denote the corresponding quantity associ¬ 
ated with electrons and ions respectively. The longitu¬ 
dinal motion is considered to be along the laser propa¬ 
gation direction (z), however transverse motion is asso¬ 
ciated along the direction of laser polarization. 
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2. Using ID approximations the above equations can be 
written as (z is direction of laser propagation). 


dz^ 


= -p/eo 


(A6) 


from Eq. A2 to A5 after separating perpendicular and 
longitudinal component we can write [laser pulse cor¬ 
responds to vector potential A (transverse) and electro¬ 
static field will correspond to (j) (longitudinal)], 

1 

^ ^ ^ - ZeuivY) (A7) 


5. The equations in dimensionless units can be written as, 

(A17) 


d^(p _ <5^ _ 


d^a d^a WepY Zn,pY 
dz^ <3f2 Yinii 


^ 9 rz y 

= -[ne -Zn; 

dt V 7g Yinii 


dtdz 


(A18) 


(A19) 


1 d^cj) 


= -HoieneV^-ZeniV) 


3. Now Consider the Lorentz force equations, 

dP 


dt 


= ^(E-t-v X B) 


(A8) 


(A9) 


dE, 

dt 


(nA-ZnA) (A20) 

V Ye Yim^ 


pY = a ; pY = -Za 


(A21) 


dP 

dt 


= q 


dA 

- V0 - — -f V X (V X A) 


(AlO) 


d^a d^a fUe Z‘^ni\ 

= ~ + — U* 
V y. Yinii) 


dz} <3f2 V Ye Yi'^^‘ 


(A22) 


dP 

dt 


= q 


r d(j) 

dA 

dA / 

[ dz 

dt 



dA\ 


(All) 


again comparing longitudinal and transverse compo¬ 
nents. 


Y^, = \+a^ + {plf-Y^ = l + ^a^+(^)^ (A23) 

mf V nii / 


dP\ 


= - — + (A24) 


/ dA 

dA\ 

dA 




dpi dtp 

/ 1 <9a2 

.Pi 


(-*- 

1 

' 

II 

; 1 

1 

II 

^P± = 

-qA (A 12) 


dt dz 

V 27 e dz 

Ye 

dz / 


—qA 

1 — 

p 2 

(A13) 

dpi 


' Z 2 da^ 

A 


v± = 

niY 

; 7=yi + 

nfic^ 

dt 

dz ^ 

. 27 m,' dz 

Yiiiii 

dz ) 


(A25) 


dt 

dP, 


= -q^+q[\^- — \ 


dz 

dtj) 


(9A\ 
dz / 


dt dz 2mY dz 


d^ dPi- 

dt my dz ^ dz 2mY dz 


dPz , Pz dPz 


(A14) 


(A15) 


(A16) 


for electrons q = —e and for ions q = Ze. 

4. We will be using following dimensionless units, 

(a) a = eAKpiec) ; (p = e^/{mec'^) 

(b) Hej = nejjnc ; Wc = £o(0^nie/c^ (critical density) 

(c) p = P/ (niec) X = kx t = at ; v = v/c ; m = 
nilnie', q = qje 


From Eq. A23 it can be proved that. 


dz 2Ye dz Ye dz 


dYi z2 da^ p\ dp[ 
dz 2Yimf dz ym? dz 

equations for pz will then be written as. 


dpi 

dt 


dYe dpi 
dz ’ dt 


dY 

' dz 


(A27) 


(A28) 


The above set of equations can be closed with continu¬ 
ity equations. 


due dr pl\ dni d , , 


dt dz 


V Ym) 




















9 


6. Now final closed set of equations to be solved are, 


Je = 




1 +3^ + (p|)2 


(A37) 


^ = Zn,- - He (A30) 

az 


7i = 




(A38) 


dt 


ne — — Zn 
^ Ye ■ Ym 


(9^a (9^a /n^ Z^ni\ 
dz? dfi We Yi^i ' 


(A31) 


7. If we consider the steady ionic background, then the 
closed set of equations reduces to, 

dE, 

= Znt - He (A39) 

az 


(A32) 


dE, 

dt 


z Pz 

= np — 


(A40) 


rz ^Ye 

dt " dz 


(A33) 


d^a d^a He 

dz^ dfi Je 


(A41) 


dt 


7TJ 

= LEr — rui -z— 

oz 


dUe 

dt 


+ 



= 0 


dnj 

dt 



(A34) 


(A35) 


(A36) 


^pI _ p ^Ye 
dt " dz 



Ye = 




1 + 3 ^ + ( p |)2 


(A42) 

(A43) 

(A44) 
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